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We analyze the stochastic evolution and dephasing of a qubit within the quantum jump (QJ) 
approach. It allows one to treat individual realizations of inelastic processes, and in this way 
it provides solutions, for instance, to problems in quantum thermodynamics and distributions in 
statistical mechanics. As a solvable example, we study a qubit in the weak dissipation limit, and 
demonstrate that dephasing and relaxation render the Jarzynski and Crooks fluctuation relations 
(FRs) of non-equilibrium thermodynamics intact. On the contrary, the standard two-measurement 
protocol, taking into account only the fluctuations of the internal energy U, leads to deviations in FRs 
under the same conditions. We relate the average (e _,3C/ ) (where /3 is the inverse temperature) with 
the qubit’s relaxation and dephasing rates, and discuss this relationship for different mechanisms of 
decoherence. 


Decoherence of quantum systems is a topic of consider¬ 
able fundamental interest, but it also presents an impor¬ 
tant technological challenge in the pursuit of quantum 
computing based on two-level systems (TLSs), qubits, 
see, e. g., p] and references therein. Quite generally de¬ 
coherence is considered on the level of averages over many 
experiments: in this situation results predicting expecta¬ 
tion values can be obtained theoretically, for instance, us¬ 
ing master equations for the density matrix of the system. 
It is, however, of considerable interest to obtain direct ac¬ 
cess to decoherence processes from the point of view of 
non-equilibrium statistical mechanics and thermodynam¬ 
ics for the system and its environment. This problem can 
be tackled conveniently by various approaches of stochas¬ 
tic quantum mechanics, e.g., via the analysis of quantum 
jumps (QJ) and non-Hermitian Schrodinger equation [2). 
It is the purpose of this article to include dephasing in 
such analysis, and to discuss various models of decoher¬ 
ence from the thermodynamic point of view. 

Work W done on a system by a source splits into the 
“useful” work, which is the change U in the internal en¬ 
ergy of the system itself, and into heat dissipated into 
the environment, Q , such that 

W = U + Q. (1) 

In order to relate such a process, for instance, to the 
thermodynamic fluctuation relations [3H5] , see [3 [7] for 
a review, the full work ( W ) needs generally to be con¬ 
sidered. There are various ways to measure work-related 
quantities in a quantum set-up. The original proposal 
was put forward as the so-called two-measurement pro¬ 
tocol (TMP), where the state of the system is measured 
first before the work is applied, and second after the ap¬ 
plication of this work BB - This yields naturally the 
difference in the internal energy in form U = Ef — E^, 
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FIG. 1. (a) The pulses applied to the qubit in TMP protocol. 

The qubit state is measured in the beginning and at the end. 
In between, two 7t/2 pulses around *-axis are applied. The 
three time intervals between the pulses are denoted n, T 2 and 
T 3 . The corresponding evolution of the qubit is shown on the 
Bloch sphere for two consequent 7r/2-pulses (t 2 = 0) in panel 
(b), and in the presence of dephasing during the time interval 
T 2 in panel (c). 


where Ei,Ef refer to the energies of the states of the 
system observed in the initial and final measurements, 
respectively. For a closed system, not interacting with 
the environment during the driving period {Q = 0), this 
yields then the whole work according to Eq. (JT|) . The 
true interesting case is, however, that where the system 
is open (Q ^ 0 in general), which cannot be captured by 
the simple TMP. Because TMP, being an obvious first 
choice, has nevertheless become a common measure of 
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work in actual experiments on quantum systems m , see 
also it is interesting to see quantitatively how good 
or bad an estimate it yields for the full work in realistic 
open set-ups. To perform such an analysis is one of the 
goals of this article, and in quantitative terms we present, 
based on the QJ method p2|, results on fluctuation re¬ 
lations on U for a generic two-level system with definite 
environmental relaxation/excitation rates and dephasing 
at finite temperatures. In the case of weak dissipation, 
we obtain leading order corrections analytically. 

To obtain the governing equation, we write the evo¬ 
lution of the wave-function of a two-level system in the 
scenario which accounts for driving, relaxation and de¬ 
phasing. In the absence of jumps the wave-function obeys 
a Schrodinger like equation with a non-Hermitian Hamil¬ 
tonian H such that [2] 

m +At)) = (i - P )- 1/2 (i - ^ W) i m), (2) 

where | = a(t)\g ) + b(t) |e), p = (T t |a| 2 + rj.|6| 2 )At 
for a small time interval of observation, At. Here, T|,T^ 
are the relaxation and excitation rates, respectively, due 
to the coupling to the bath, and they obey the detailed 
balance condition T| = e~^ hu °T^ at inverse tempera¬ 
ture /3 for level spacing htoo of the qubit. These rates 
determine the probabilities of interrupting the evolution 
of Eq. ([2| by excitation, p^ = T||a| 2 At, or by relax¬ 
ation, p± = T^|6| 2 At, during the time interval At. | g) 
and |e) denote the qubit ground and excited state, re¬ 
spectively. In the interaction picture we take the unper¬ 
turbed Hamiltonian Hq of a TLS with the level spacing 
equal to the drive frequency too- Thus, we may write the 
total non-Hermitian Hamiltonian (excluding that of the 
bath itself) of the quantum trajectory approach as 

H = H 0 + SH + V + N. (3) 

Here H 0 = (fiwo/2)(|e)(e| — \g)(g\)- The slow fluctuation 
of the level spacing responsible for dephasing reads SH = 
(hSuj/ 2)(|e)(e| — |g)(g|). We assume that Su(t) <C too is 
a classical stochastic process. The driving term, with 
the drive signal A (t) is chosen to be of the form V = 
A(t)(|g)(e| + |e)(g|), and the non-Hermitian part of the 
Hamiltonian reads N = — (jfi/2)(P||e)(e| +r^|g)(g|). We 
write the Schrodinger equation for the wave-function in 
the interaction picture, |= e lHat l h \ip(t)). Then 
expanding up to linear order in At, we find the evolution 
corresponding to Eq. <§: 

| ipi(t + At)) = (1 - p)~ 1/2 x 

l-UVj^+SHj^ + N^At 

h 


I Mt)h (4) 


tion of a and b of = a\g) + b\e) obeys then 

a = —^A (t)e~ lUot b + i^a + ^-a|&| 2 , 

b=- lA(t)e- ot a-^6-^|a| 2 &. (5) 

Here we have defined AT = Tj, — T-f. 

In what follows, we analyze an exemplary protocol in 
the weak dissipation regime, and demonstrate that, as 
expected, the Jarzynski equality (JE) 

(e~ pw ) = 1 (6) 


holds, whereas (e~^ u ) differs from unity. Before present¬ 
ing the formal derivation of the weak dissipation results, 
we justify such a deviation by a simple argument in the 
limit where the time interval between the two measure¬ 
ments is very long. In general, we can write the TMP 
outcome as 


=pU 4 
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(7) 


Here p l a is the probability of detecting state a = g,e at the 
first measurement and p^ a is the conditional probability 
of measuring state b at the end (2nd measurement, /), 
if the system was initially in state a. Taking the time 
interval between the two measurements to be very long 
means that all the populations in Eq. Q are thermally 
distributed, p’ = 1 - p\ = p f g]g e = 1 - p f e[g e = (1 + 
e -phu, 0 )-i 5 anc i Eq _ 0 yields 

(e ^ l3U 'j = 1 + tanh 2 (/Ifia;o/2). (8) 


In the high-temperature limit (/ 3Hu>o <C 1), this ex¬ 
pression yields (e~^ u ) = 1, but at low temperatures 
(/3fuoo 1), we obtain {e~ tiu ) = 2 in stark contradic¬ 
tion with JE. This result holds for any driving protocol 
between the two measurements. 

Now we move to the weak dissipation treatment tak¬ 
ing into account trajectories which include at most one 
relaxation or excitation event. This way we obtain devi¬ 
ations from the fluctuation relations linear in the transi¬ 
tion rates. We may write for the work exponent 

(e~ fiw ) = Po(e- /3w )o + Pi{e~ 0w )i, (9) 


where Pq, Pi are the probabilities of zero and one-photon 
processes, and (e~P w )o, (e~^ w )i are the averages for the 
corresponding processes. We will explicitly check the va¬ 
lidity of JE up to this order. On the other hand, we have 
for the quantity measured in the TMP protocol, 

(e~P u ) = P 0 (e-P u )o + Pi(e-^)i (10) 


where operator O in the interaction picture reads Oi(t) = with the corresponding notations. Now for zero-photon 

e iH a t/hQ e -iH 0 t/h _ Eor the two state system, the evolu- processes, U = W, thus (e -/3c/ )o = ( e~^ w )o■ 
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We choose the following protocol, see Fig.[T| The qubit 
is measured in the beginning and at the end of the pro¬ 
cess (TMP), and it is driven in between by two consec¬ 
utive pulses, each of them producing ideally a 7r/2 rota¬ 
tion around the same, say x-axis. We assume that the 
pulses are short enough in time such that no relaxation 
or dephasing takes place during them. This is possible to 
achieve by adjusting the drive amplitude of the rotation 
pulses. The timing outside the pulses is such that the 
first 7 t/ 2 pulse is applied at time T\ after the initial mea¬ 
surement of the qubit state. The time interval between 
the two rotation pulses is denoted 72 , and during this 
time dephasing can take place. Finally, T 3 is the time in¬ 
terval between the second rotation pulse and the second 
measurement pulse. This sequence has the property that 
without relaxation and dephasing it performs a standard 
7 r rotation swapping the system between the states g and 
e. Yet, the qubit may jump at any moment during the 
protocol, if F’s are finite, and in between the two tt/2 
pulses the qubit is in a superposition state and it is sus¬ 
ceptible to dephasing as well. Typical trajectories on the 
Bloch sphere are illustrated in Fig. [l] (b) for r 2 = 0 ( 7 r 
pulse) and in (c) for r 2 > 0 (dephasing) starting in the 
state g at the north pole. 


Derivation of the averaged work exponents (18 1 and 


(19) is made along the following lines. We follow the 


amplitudes along the operation of the qubit according to 
the protocol shown in Fig. [T| We take into account those 
trajectories that involve either no or one quantum jump. 
Knowing the amplitudes we evaluate the probabilities of 
different trajectories and then calculate the thermal av¬ 
erages up to the first order in T’s. The detailed cal¬ 
culation, outlined in the Supplemental Material, Sec. 1, 
results in the relationship (e~^ w )s u = 1 where the sub¬ 
script Suj emphasizes that we average over all realizations 
with a fixed level spacing ui 0 + Suj. This result implies 
that the Jarzynski equality is valid for any distribution 
of Suj, {e~P w ) = 1 . On the contrary, for U we obtain 

(e~ /3U )s ul -l = [r 3 -TiCOs(^ 2 )]r s tanh 2 (/3Swo/2) (11) 

where F E = F^ + F| and Sp 2 = f^ +T2 Suj(t) dt. After 
averaging over realizations of the random process Su>(t) 
we get 

(e~' 3U )-l = [T3-ri(cos(^ 2 ))]r s tanh 2 (/3?kuo/2). (12) 


Equation (12) is the central result of the present work 
it relates the difference A = (e~P u ) — 1 with decoherence 
of the qubit characterized by the average (cos (S(p 2 ))- If 
the qubit is fully characterized and the above average is 
known, then Eq. ( |l2| ) can be used for checking the JE. 
Indeed, since (e - ^” 7 ) = 1 we can correct for usage of U 
instead of W. If r E (ri +t 2 +t 3 ) <C 1, then the correction 
is negligible. 

On the other hand, if the quantity A({ri}) is deter¬ 
mined from the experiment (say, by direct registration 


of emitted or absorbed photons), then an information 
regarding dephasing can we extracted as 


(cos(<5</> 2 )) 


73 1 _ / din A 

7i \ d In 73 


-i' 


(13) 


The quantity (cos(<5 ip 2 )) contains the same information 
as what can be obtained in a Ramsey measurement |13j . 
This quantity appears due to the specially selected pro¬ 
tocol shown in Fig. [I] Other protocols would contain 
different averages, and in this way may provide an ad¬ 
ditional information about the underlying decoherence 
mechanism. For instance, when a 7r/4-rotation is first 
performed around ?/-axis, and then a 7r/2-rotation around 
x-axis, one obtains corrections oc (sin(5</j 2 )). 

In the absence of dephasing (cos(<fy? 2 )) = 1 and 

(e~ pu ) - 1 = ( 7-3 - ri)r E tanh 2 (/3Sa;o/2). (14) 

The result depends on (73 — iq) through the difference of 
the probabilities of photon-assisted processes during the 
time slots 73 and t -\ . respectively. 

Figure [2] shows the analytic predictions and numer¬ 
ical simulations of the quantities analyzed above for a 
chosen set of parameters as detailed in the figure cap¬ 
tion. The offset Suj is assumed to be constant for each 
set of data. The numerical calculations are performed 
using a stochastic simulation with the QJ method, as¬ 
suming ideal rotations. What one observes apart from 
the statistical scatter due to the finite number of repe¬ 
titions in the averages, 10' for each point, is that the 
analytic approximations follow closely the numerical re¬ 
sults. Only for long r 3 intervals in (a), we see that the 
linear approximation for (e~P u )s u overestimates the de¬ 
viation from unity. This is natural, as we have shown 
that for 73 —► 00 , (e~ f3U )sui ~ 1-97 for the parameters in 
Fig. [ 2 ] see Eq. (|8|. On the contrary, in all cases (e~^ w )s u 
stays around unity within the statistical scatter. 

One obtains further insight of the differences between 
U and W by inspecting the actual probability distribu¬ 
tions of internal energy and work under the given driving 
protocol. Both ratios U/Hujq and W/Hujq can obtain in¬ 
teger values. Figure [5] shows the numerically calculated 
distributions for U and W under the same conditions as 
in Fig. [2] varying the delay time r 3 between the second 
rotation and the final measurement. There are several 
important conclusions to draw from the dependencies in 
Fig. i Firstly, the probabilities pu(U) relax with time 
73 since U measures the internal energy which changes 
after the non-equilibrium driving. On the contrary, the 
waiting period 73 does not influence pw( W) since there 
is no work done after the second rotation pulse. These 
dependences are consistent with our basic expectations. 
Finally, we also recover the important Crooks fluctuation 
relation for work [5] in the form 

Pw(W)/ Pw (-W)=ef )w (15) 
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for those values of work that are within the reach in this 
situation (for W/hojQ = ±1). On the contrary, the ratio 
Pu(U)/pu(—U) is not constant in 73 and thus does not 
satisfy a fluctuation relation. 

Neglecting time dependence of 6u> during the time in¬ 
terval T 2 and assuming a Gaussian distribution of 5u> with 
vanishing mean and width (Soj 2 ) = 2T^ p we obtain 

(e~ /3C/ ) — 1 = [(73 — e~( r ' pT2 ^ 2 ti] Ts tanh 2 (/3?kuo/2). (16) 

Therefore, A is directly related to T^, which is the decay 
rate for (cos(5^2))- This is the dephasing rate of the 
qubit during the time interval 72 . 

A realistic source of decoherence caused by a low- 
frequency noise is discussed in the Supplemental Mate¬ 
rial, Sec. 2. We consider so-called random telegraph noise 
created by a two-state dynamic degree of freedom in the 
environment, as well as l//-noise produced by overlap 
of many degrees of freedom, see mm for a general re¬ 
view and [6] for a review on qubits. Experiments [5] 




FIG. 2. Analytic expressions given in the text, compared 
to numerical results obtained by stochastic simulations, see 
Ref. [12] for details. The protocol is the one described in the 
text and in Fig. IT] We assume /3huo = 5 and F^/Awo = 0.03 
(circles), and r^/ftwo = 0.1 (triangles). In (a) won = W 0 T 2 = 
1, wo 7"3 varies, and Su> = 0. The open symbols refer to 
( e -/ 3 W)^, an£ j t j le gjj ec j oneg to (e _flc/ ) 5lJ . The solid lines 
represent the analytic predictions of Eq. (Ill for {e~ l3U }su> 
for the two values of Fj.. In (b) {e~ l3U )su> is plotted for 
won = won = W 073 = 1 against Sifi 2 = 5lut2 with other pa¬ 
rameters as in (a) (solid lines and filled symbols). The results 
for ( e~ l3W )sui , indicated by the corresponding open symbols, 
are again concentrated around unity. In all cases we employed 
10 ' repetitions of the protocol for each data point. 



FIG. 3. Probabilities of different possible outcomes of u = 
U/huJo and w = W/Hujq. We assume the same protocol as 
before with /3hu>o = 5, Ti/hujo = 0.1, 5ui = 0, and won = 
wo 72 = 1. We vary the delay time of the second measurement, 
73. The solid lines refer to U and the dashed ones to W. We 
employed 3-10 5 repetitions of the protocol for each data point. 


on Josephson qubits have indicated that charged defects 
may be responsible for the 1// noise. 

Another mechanism of decoherence is a back-action of 
the measuring device, see, e.g., m- In this case the de¬ 
coherence rate can be related to the so-called accessible 
information , i.e., the information gained in the measure¬ 
ment process, see, e.g., [T5] and references therein. If this 
decoherence mechanism is dominant, Eq. (12) allows re¬ 
lating thermodynamics of the qubit to the accessible in¬ 
formation gained by the detector. Note, however, that 
the exact relationship between the decoherence and the 
accessible information depends on the detector type and 
the measurement protocol. 

The decoherence and thermodynamics can be related 
to the mutual information also in the case when decoher¬ 
ence is due to interaction with fluctuators, but only when 
one is able to measure both the qubit and the fluctuators. 
The relationship between decoherence and mutual infor¬ 
mation has been analyzed for a qubit interacting with a 
single fluctuator coupled to phonon environment |19| . 


To summarize, we have studied the distribution of both 
the internal energy U and the work W done on a qubit 
in a driven non-equilibrium protocol under weak dissipa¬ 
tion. Specifically we relate the non-equilibrium thermo¬ 
dynamic relations to decoherence rates. We demonstrate 
that the common fluctuation relations are satisfied for 
W, but not for U which is the quantity measured in the 
two-measurement protocol, TMP. Our results allow one 
to evaluate the validity of, e.g., the Jarzynski equality 
from measurements of the internal energy U rather than 
of the total work W. From the opposite point of view, if 
the average {e~^ u ) is measured, one can use the results 
for studying dephasing of the qubit via such a statistical 
measurement. 
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Dephasing and dissipation in qubit thermodynamics: 
Supplemental material 

1. Jarzynski equality in the weak dissipation limit 

Hereby we derive the expression for (e~P u ), Eq. (11) 
of the main text, and prove that the Jarzynski equality 



tributing to various Pi{-)i averages described in the text. On 
the left of each diagram, we indicate the result of the ini¬ 
tial measurement ( g, e), the vertical dashed lines indicate the 
rotations of the qubit and on the right we indicate the out¬ 
come of the final measurement ( g,e ). The relaxation and 
excitation events are indicated by down and up arrows, re¬ 
spectively. The jump occurs in the time interval 0 < t < i~i 
in (a), n < t < n + t 2 in (b), and n + t 2 < t < n + r 2 + t 3 
in (c). 


(JE) 


{e ~P w ) = i 


(17) 


holds in the weak dissipation regime. In this limit we 
take into account trajectories which include at most one 
relaxation or excitation event, see Fig. [4j and obtain cor¬ 
rections to fluctuation relations linear in the transition 
rates. We may write the work exponent as 

(e-? w ) = P 0 (e~ fiw )o + Pi(e~P w )i, (18) 


where Pq, Pi are the probabilities of zero and one-photon 
processes, and (e~^ w )o, (e~^ w )i are the averages for the 
corresponding processes. We will explicitly check the va¬ 
lidity of JE up to this order. On the other hand, we have 
for the quantity U measured in the two-measurement pro¬ 
tocol (TMP), 

(e-W) = Po(e~ pu )o + Pi(e~ pu )i (19) 


with the corresponding notations. Now for zero-photon 
processes, U = W, thus (e~P u ) 0 = (e~^ w ) 0 . Following 
the protocol outlined in Fig. 1 of the main text we find 
the amplitudes along the operation trajectories of the 
qubit. If the system is in state g just before the first it/2 
rotation, it acquires amplitudes 

a gi T i+) = 1/ V%, b g (r 1+ ) =l/V2 

right after it. Similarly 

a e (ri + ) = — 1 /a/2, b e {Ti+) = l/\/2. 


In the absence of jumps, the system evolves during the 
time interval r 2 such that according to Eq. (5) of the 
main text, up to linear order in AT, 


a g (t 2 -) = -a e {t 2 -) = -^=(1 + Arr 2 /4)e I ' 5 ‘ P2/ ' 2 , 
b g (t 2 .) = b e (t 2 -) = -^(1 - AIY 2 /4) e -^/ 2 , 
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where t 2 = t± + 72 , Sip 2 = / Q r2 Soj(t) dt. After the second 
7 t/ 2 rotation in this situation, 

a g (t 2+ ) = -b e {t 2 +) = i sin(( 5 v 3 2 / 2 ) + ^ Arr 2 cos(6tp 2 /2), 

6g(t 2+ ) = —a e (t 2 +) = cos(<5<p 2 /2) + - ATt 2 sin((j</? 2 /2). 

Finally, the probabilities at the moment of the second 
measurement, £3 = T\ + t 2 +T 3 , again assuming no jumps 
after the first rotation, are given by 

I a g (t 3 )\ 2 = sin 2 (6(p 2 /2) [1 + Arr 3 cos 2 (^ 2 / 2 )], 
|a e (f 3 )| 2 = cos 2 (<fy 2 /2)[l + Arr 3 sin 2 ((5v? 2 /2)], 

with | 6 Si e(t 3 )| 2 = 1 - \a g , e {t 3 )\ 2 - 

With the help of the evolution of the amplitudes, we 
can construct the averages. In particular for no-jump 
trajectories we may write 

P 0 (e~P u )o = Po(e~P w )o = 

Pg e- d ^\\a g (t 3 )\ 2 e° + \b g (t 3 )\ 2 e-^) 

+p e e~ d ^\\a e (t 3 ) + |& e (t 3 )| 2 e°). ( 20 ) 

Here p g = 1 — p e = (1 + e - ^" 0 ) -1 , 
di(r)= 

Jo 

Taking only contributions yielding up to linear correc¬ 
tions in T’s, we obtain 

Po(e- 0u )o = Po(e- /5w )o = l- 

p g T i Ti[cos 2 (Sp 2 /2){l + e' 2l3hu °) + 2sm 2 (Sp 2 /2)e~ /3hu10 } 
-\^T 2 -2p g T i T 3 e-^. (21) 

Here Ts = + T-p For the single-jump trajectories, 

we follow the diagrams in Fig. [4] in this supplemental 
material to obtain again the contributions up to linear 
order in F’s. As an example, for the diagrams in Fig.[4](a) 
we obtain 

Pi,a(e~ PU )\,a =P 9 r t Ti(|a e (t 3 )| 2 e 0 + \b e {t 3 )\ 2 e~ phu °) 

+Pet'iT 1 (\a g (t 3 )\ 2 e phu>0 + \b g (t 3 )\ 2 e°) 

= PgF^n[sin 2 (<5^/2) (1 + e" 2 ^°) ( 22 ) 

+2 cos 2 (5<p 2 /2)e~^ hu °], 


and 

Pi, a (e~ 0W )i,a = p g T tTl (\a e (t 3 )\ 2 e^ + \b e (t 3 )\ 2 e°) 

+p e r i ri(|ag(t 3 )| 2 e° + \b g (t 3 )\ 2 e~ phui0 ) 

= p g Y±T\ [cos 2 (5<p 2 /2)(l + e ~ 2 P hu °) (23) 

+2sin 2 ((5^ 2 /2)e- /3 ^ 0 ]. 

In these equations Py^ri « 1 —e _r ^' tTl is the probability 
of a jump within the time interval t\. 


For the diagrams in Fig. [4] (b) we obtain 

Pi,b{e~ 0u )i :b = 

Pg I 't\a g ,-K/2\ 2 r 2 {\a et „ /2 \ 2 e 0 + \b e ^ /2 \ 2 e~ 0hu °) 

PPg^i\^ g ,Tv/2\ 2T 2{\cL gi7r / 2 \ 2 e° + |&g j7r /2| 2 e 

+p e r-|-|a eiW / 2 | 2 r 2 (|a eiW / 2 | 2 e /3?icJo + \b et7r / 2 \ 2 e°) 

+Pe^'l\b e ,TT/ 2 \ 2 ' T 2 (\ 0 ’ gt T T / 2 \ 2 e^ hUo + |£>g i7r /2 | 2e °)- (24) 

Here \a g ^ /2 \ 2 = \a e ^ /2 \ 2 = \b g ^/ 2 | 2 = \b e ^/ 2 \ 2 = 1/2 are 
the probabilities after a n/2 pulse when the system starts 
from an eigenstate. Inserting these values we obtain 

Pi,b(e~ PU )i,b = ^eT2. (25) 


Similarly, 

Pi,b(e~ 0W )i,b = 

Pg T t\ a g,ir/2\ 2 T 2 (\a etn / 2 \ 2 e 0huo + |6 e>7 r/ 2 | 2 e°) 

PPgFl\b g ,7r/ 2 | 2 ' r 2 (|ag,7r/ 2 | 2 e + \b g ^/ 2 \ 2 e 2 P hu °') 

+PeHa e ^ /2 \ 2 T 2 (\a e ^ /2 \ 2 e 2 ^ + \b e , n/2 \ 2 e^°) 
+P e r4.|&e,7r/2| 2 ' r 2(|ag,7r/2| 2 e° + |^g,7r/2| 2 e ,3fi ' aJo ), (26) 
yielding 

Pi,b(e~ 0W )i,b = ^r E r 2 . (27) 

For the diagrams in Fig. [4] (c), we have 
Pi,d e ~ PU )i,c = 

p g (T t \a g (t 2 + )\ 2 T 3 e-^° + r 4 .|6g(t 2+ )| 2 r 3 e 0 ) 
+p e (r t |a e (t 2+ )| 2 r 3 e 0 + r ; |6 e (t 2+ )| 2 r 3 e /3fiaJo ), (28) 

yielding 

Pi,c( e ~ PU )hc = Ps^M 1 + e- 2 ^°). (29) 

Similarly we obtain for W for the diagrams in Fig. [4] (c) 
Pl,c( e ~ PW )l,c = 

Pg(r t |ag(t 2+ )| 2 r 3 e 0 +r ; |6g(t 2+ )| 2 r 3 e-^) 
+Pe(J'^\a e (t 2 +)\ 2 T 3 e /3hu ° +r 4 .|& e (t 2+ )| 2 r 3 e 0 ), (30) 
which gives 


Pi,c(^ pW )i,c = 2p g T i T 3 e-^_ 


(31) 


Combining Eqs. (21), (23), (27) and (31), we obtain 
up to the first order in T’s, 


(e-P w ) Su = Po(e~ f3w )o+ J2 p iAe- 0W )i,i = 1, (32) 


i=a,b,i 


where the subscript dui emphasizes that we average over 
all realizations with a fixed level spacing uiq + Sto. This 
result implies that the Jarzynski equality, (e~J w ) = 1, 
is valid for any distribution of 6to. Similarly, combining 


Eqs. (21), (22), (25) and (29), we obtain Eq. (11) of the 


main text. 
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2. Decoherence due to low-frequency noise 

Here we will discuss a realistic source of decoherence 
caused by a low-frequency noise. A part of these fluctu¬ 
ations typically has a 1 // spectrum and is referred to as 
1 // noise. Such noise is generic for all disordered mate¬ 
rials (for a review see, e. g., [I] and references therein). It 
is also common in single-electron and other tunneling de¬ 
vices, see, e. g., [2]. Experiments [3] on Josephson qubits 
have indicated that charged impurities may be responsi¬ 
ble for the 1 // noise. 

One of the most common sources of low-frequency 
noise is the rearrangement of dynamic two-state defects, 
fluctuators , see, e. g., [3] and references therein. Ran¬ 
dom switching of a fluctuator between its two metastable 
states (1 and 2) produces random telegraph noise. The 
process is characterized by the switching rates 712 and 
721 for the transitions 1 — > 2 and 2 —> 1. Only the fluc¬ 
tuators with energy splitting E < T, contribute to the 
dephasing since the fluctuators with large level splitting 
are frozen in their ground states (we measure tempera¬ 
ture in the energy units). As long as E < T the rates 712 
and 721 are close in magnitude, and without loss of gen¬ 
erality one can assume that y 12 « 721 = 7 . i. e., the fluc¬ 
tuations can be described as a random telegraph process , 
for reviews see mm- A set of random telegraph fluc¬ 
tuators with exponentially broad distribution of relax¬ 
ation rates, 7 , produces noise with 1 // power spectrum 
at 7min <C oj = 27r/ < 7 0 . Here 7 m j n is the switching 
rate of the “slowest” fluctuator while 70 is the maximal 
switching rate for E = T. Random telegraph noise has 
been observed in numerous nanodevices based both on 
semiconductors, normal metals, and superconductors 0 . 

For evaluating role of fluctuators in dephasing, we will 
use a simple classical model within which one can ana¬ 
lyze exactly the qubit response to typical manipulation 
protocols. According to this model, the quantum system 
- qubit - interacts with a set of two-level entities, see [ 6 j 
for a review. The latter fluctuate stochastically between 
their states due to interaction with a thermal bath, which 
may be not directly coupled to the qubit. 

We start from the essentially non-Gaussian situation 
when there is a single fluctuator coupled to the qubit. 
Following [B], we assume that the fluctuator is described 
by the Hamiltonian 7-Lp = {E/2)t z where t z is the Pauli 
matrix of the fluctuator and E = VA 2 + A 2 is its energy 
splitting. The latter depends on the diagonal splitting A 
between the fluctuator’s states and their tunneling cou¬ 
pling A. 

This fluctuator switches between its states due to the 
interaction with environmental bosons (phonons or elec¬ 
trical fluctuations), the switching rate being [Bj 

7 = (A/E) 2 7 o(E). (33) 

A typical estimate for 70 is T 3 /HT 2 , where T* ~ 20 K. 


Following the approach outlined in [ 6 ], we represent 
5w(t) as a random telegraph process, x(t), as Sw(t) = 
v\(t). Here x(f) switches between the values ±1 at ran¬ 
dom times, distributed according to the Poisson distri¬ 
bution, (x(t)x(O)) = e~ 2lt . The coupling constant v, 
in general, depends on the operation point of the qubit, 
v oc (A q /Huio)(A/E), so the mechanism can be verified 
by changing this point. Here A q is the tunneling splitting 
of the qubit. 

The average 7 |r 2 ) = (cos(<5 ip 2 )) reads as 0 

p--yr 2 

[(m + l)e WT2 + {p- l)e _7MT2 ] (34) 

where p = \J\ — (u/ 7 ) 2 . At v -C 7 , 

i>(v, 7 N) ~ e~ r * T2 , Tq, = v 2 /2-y. (35) 

This behavior is typical for the so-called motional nar¬ 
rowing of the spectral lines in magnetic resonance [7j. At 
t)>7, 


V’( v i7l r 2) ~ e 7T2 cos('CT 2 ) , (36) 


i.e., we have beating between the qubit’s levels Eq ± v. 
These dependences are shown in Fig. [5} 

Now we consider the case of many fluctuators produc¬ 
ing 1 // noise. We will use the same approach as in [ 6 ] 
and assume that the random processes of different fluc¬ 
tuators are not correlated and their total number N 3> 1. 
Then (cos(< 5 <^ 2 )) can be expressed as e~ K '( T2 ' > where 


K-{t 2 ) = J dvd'yV(v,'y) [1 - ip{v,flT 2 )\ ■ (37) 


Here <ij}(y,^\t) is given by Eq. (34), while V(v,j) is the 


distribution of the coupling parameters, v, and the fluc¬ 
tuator’s switching rates, 7 . 

To specify the distribution function one has to assume 
a model for fluctuators. An overview of different models 
is given in Refs. [6j. To formulate the results let us in¬ 
troduce a typical coupling strength rj as the interaction 
strength at the average distance between fluctuators with 
the interlevel spacing of T. Here we will briefly discuss 
the case when fluctuators are uniformly distributed in 



FIG. 5. Plots of the function ip(v,'y\i) given by Eq. (341 for 
v/y = 0.5 (blue curve) and v /7 = 5 (brick-red curve) . 









the space around qubit and the interaction between the 
qubit and a fluctuator decays <x 1/r 3 . In this case [ 8 , 9] 


/C(r 2 ) « rj ■ 


7 o t 2 for 7oT 2 < 1; 

r 2 ln 7 0 r 2 for 7 0 r 2 1 . 


(38) 


The results for the first limiting case coincide with the 
Gaussian decay, Eq. (16) of the main text, with « 
y/fp/Q. This result has a clear physical meaning: the de¬ 
coherence occurs only provided that at least one of the 
fluctuators flips. Each flip provides a contribution ~ ? 7 t 2 
to the phase, while 70 r 2 <C 1 is a probability for a flip 
during the observation time. The result for 70 t 2 1 is 
less intuitive since in this region the dephasing is non- 
Markovian, see [5j for more details. 

It is important that at large observation times, r 2 
7 ( 7 1 , the decoherence is dominated by few optimal fluc¬ 
tuators. The distance r opt (T), between the optimal fluc¬ 
tuators and the qubit is determined by the condition 


v{r opt ) « 70 (T). 


(39) 


This estimate emerges naturally from the behavior of the 
decoherence in the limiting cases v 7 and 1 K 7 . For 
strong coupling the fluctuators are slow and the qubit’s 
behavior is determined by quantum beatings between the 
states with E ± v. Accordingly, the decoherence rate is 
of he order of 7 . In the opposite case, as we already 
discussed, the decoherence rate is oc v 2 / 2y. Matching 
these two limiting cases one arrives at the estimate (391. 


Let us discuss possible implications of the results for 
the suggested TMP. We start with the case of a single 
strong fluctuator. In principle, one can tune the qubit 
(i.e., the manipulation frequency) to (Soj) = 0 using 
the TMP. Note, however, that such a tuning would be 
temperature-dependent since (Su>) generally depends on 
temperature through the fluctuator’s occupation num¬ 
bers. Therefore, it may happen that tuning should be 
made for each temperature. 

As we have discussed, the behavior of ^(r, 7 |r 2 ) de¬ 
pends on one dimensionless parameter, u /70 (T). In both 
cases theory predicts exponential decay, that corresponds 
to the Lorentzian spectrum. However, the decay rates in 
the different regimes are very different depending on tem¬ 
perature. In the case of a “weak” fluctuator, the decay 


rate is P v = v 2 /2 7 , cf. with Eq. (35), where 7 oc 70 (T), 


see Eq. (33), is increasing with temperature (the typical 


dependence is oc T 3 ). 

At low temperatures one can expect the parameter 
v/"/q(T) to become large. In this regime one can expect 
damped beatings with frequency v, see Eq. (36), the de¬ 


cay rate being 7 . In principle, one can find the coupling 
strength, v, from the beatings’ frequency. If a crossover 
between the regimes is experimentally feasible, then one 
can be sure that at the crossover point the fluctuator 


is optimal in the sense of Eq. (39), i.e., it is the most 
harmful. 


If the decoherence is produced by many fluctuators 
(l//-noise) we have to study Eq. (38), which also pre¬ 
dicts a temperature-dependent crossover at 70 (T)t 2 ss 1. 
At low temperatures and relatively small delay times, 
70 (T)r 2 <C 1, the predicted decay is Gaussian. The pre¬ 
dicted decay rate in this case, T^, « y/rj{T)^ 0 {T) oc T 2 . 
Note that rj is proportional the number of fluctuators 
with E < T, therefore rj{T) oc T. At higher tempera¬ 
tures a crossover to the case 7o(T)t 2 1 takes place. 

In this - non-Markovian - regime, the predicted decay is 
exponential, the decay rate being ij{T) oc T. 


The above discussion holds for the situation when the 
fluctuators are evenly distributed in a 3D device and the 
interaction between the qubit and the fluctuators decays 
oc 1/r 3 . If the fluctuators are formed close to some low- 
dimensional surface, then the situation can be more di¬ 
verse. Typically, in that case a small group of fluctuators 
becomes dominant, the decoherence is non-Gaussian, and 
pronounced mesoscopic fluctuations of the decoherence 
time can be expected. Generally, the behavior of the de¬ 
cay can be accounted assuming that there is one or few 
strong fluctuators and applying expressions for a single 
fluctuator, see more details in [Hj. 
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